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The galactic center excess of gamma ray photons can be naturally explained by light Majorana 
fermions in combination with a pseudoscalar mediator. The NMSSM provides exactly these ingredi¬ 
ents. We show that for neutralinos with a significant singlino component the galactic center excess 
can be linked to invisible decays of the Standard-Model-like Higgs at the LHC. We find predictions 
for invisible Higgs branching ratios in excess of 50 percent, easily accessible at the LHC. Constrain¬ 
ing the NMSSM through GUT-scale boundary conditions only slightly affects this expectation. Our 
results complement earlier NMSSM studies of the galactic center excess, which link it to heavy Higgs 
searches at the LHC. 
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I. INTRODUCTION AND BASICS 


While the existence of cold dark matter as the main matter component of today’s Universe is generally 
acknowledged, the particle nature of it is still elusive. Searches for dark matter coupled to Standard Model 
fields with more than a gravitational interaction strength follow three distinct strategies: direct detection, 
indirect detection, and production at colliders. The latter will receive a significant boost with the start of LHC 
Run II. The key question is how in the case of weakly interacting dark matter the different search strategies can 
support and inspire each other. One of the main search strategies for dark matter at the LHC are invisible Higgs 
decays, most notably in weak boson fusion mm- For example, in models without new strongly interacting 
particles such invisible Higgs decays will drive mono-jet searches and are likely to dominate over dark matter 
pair production in weak boson fusion [3]. In this paper we will reinforce the link between the specific Fermi 
galactic center excess [4]-[6] and invisible Higgs decays at the LHC [71E] in the framework of the next-to-minimal 
supersymmetric Standard Model (NMSSM) [9lITO]. 

The Fermi gamma ray space telescope searches for dark matter signals in its photon spectrum. An excess of 
gamma rays from the galactic center has for many years avoided possible background interpretations. It can 
be explained by annihilating dark matter with a spherical distribution around the center of our galaxy. Its 
spectrum gives preferred mass values for different dark matter candidates. For an annihilation to bottom quarks 
the preferred mass of the dark matter agent (Hooperon) ranges around 40 GeV [11[TT], extending all the way to 
70 GeV [12]. In our analysis of the bb case we will follow Ref. [12] and assume a conservative LSP of 30 GeV to 
70 GeV. The cross section should be in the range of av ^ 1.8 • 10“^^cm^ /s ana, with appropriate theoretical 
or parametric uncertainties for example from the choice of dark matter profile, consistent with the latest Planck 
thermally averaged results [5] . Such values are intriguingly close to the expectations for a thermally produced 
weakly interacting dark matter particle (WIMP) [T3lIT^ . 


In the MSSM the preferred mass range of the Hooperon is a challenge and typically relies on dark matter 
annihilation into a pair of gauge boson na. In the absence of the highly efficient annihilation through an 
s-channel mediator decaying for example into bb pairs the predicted relic density in the Universe tends to be 
too large. Finding efficient annihilation channels is a serious issue in supersymmetric models [TBUTS] : first, 
s-channel annihilation through the Z-pole, the SM-like Higgs resonance iLi 25 , or a heavy Higgs resonance 
are either forbidden by other constraints or too small. Any co-annihilation channel requires an additional 
supersymmetric particle within 10% of the LSP mass [T9]-[2T], which is disfavored by LEP constraints [22]. The 
way out is an additional mediator, ideally a pseudoscalar with a mass not far above twice the LSP mass m- 
This feature is clearly visible in an analysis in terms of simplified models or effective field theory [24] . 


As an extension of the MSSM, the NMSSM provides exactly such a mediator, the pseudoscalar part of the 
singlet/singlino superfield mixed with the MSSM-like pseudoscalar [25]. In the required mass range it will 
naturally decay to bb pairs, and with a reduced branching ratio to r+r“. Such an NMSSM setup can be 
tested in a parameter scan [26] and then linked for example to 4-body Higgs decay [27], trilepton searches at 
the LHG [28] or even the electroweak phase transition [29]. Because of the structure of the NMSSM we can 
follow two strategies to accommodate the galactic center excess laniisi]: first, we can keep the standard bino- 
Higgsino LSP composition of the MSSM and only couple the neutralinos to the light additional pseudoscalar. 
Alternatively, we can replace the bino content by a singlino content and assume a singlino-Higgsino, or better 
bino-singlino-Higgsino LSP. Again, it will couple to the pseudoscalar mediator. 

For the bino-Higgsino case the relevant LSP and mediator states are not decoupled from the Standard 
Model. This means that for example the pseudoscalar mediator can be searched for at the LHG [30]. The 
singlino-Higgsino channel is more challenging. After introducing the NMSSM singlet/singlino sector and its 
phenomenology including some useful formulas, in Sec. [HA] we will study its TeV-scale parameter space linked 
to the galactic center excess. In Sec. HB we will link the Hooperon parameter space to the size of invisible 
Higgs couplings. It will turn out that similar to a dark matter Higgs portal [32] the NMSSM interpretation 
of the galactic center excess will lead to invisible Higgs decays with branching ratios accessible during the 
upcoming LHG run. In Sec. |HI| we will apply the same criteria to a high-scale NMSSM setup. For this model a 
global SFitter likelihood analysis is useful, before we turn to the link between the galactic center excess and 
invisible Higgs decays. 
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A. NMSSM 

Compared to the minimal supersymmetric Standard Model the superpotential of the NMSSM [3110] includes 
an additional singlet superfield S and the associated terms 


W^NMSSM = Wmssm + A SHuHd + ^fS + !^S^ + '^S\ (1) 

where A and hz are dimensionless couplings coupling the singlet to itself and to the Higgs bosons. When 
the singlet acquires a vacuum expectation value Vg, the Higgs-singlet mixing introduces an effective /i-term 
/^eff = A'T’s. The quadratic term proportional to ji' is the supersymmetric mass term for the singlet, comparable 
to the /i-term for the MSSM Higgs bosons. Assuming a global supersymmetry the tadpole term proportional to 
can be removed though a constant shift of the singlet field. Finally, with the help of an ad-hoc Zs-symmetry 
we can make the superpotential scale invariant and set the one remaining dimensionful parameter, /i', to zero. 

The extended superpotential in Eq. 0 in terms of the superfield S can be translated into additional soft- 
SUSY-breaking terms for the physical singlet field 5 m, 

5'^ + + h.c.^ . ( 2 ) 

The Ax,k carry mass dimension and fix the scale of the Lagrangian, while A and k, defined in Eq. 0 are c- 
numbers. An alternative parametrization of the same Lagrangian uses the mass terms = Bjj. and rug = H'/i'. 
To be consistent with the Zs-symmetry of the superpotential we also eliminate the corresponding SUSY- 
breaking terms by setting m 3 = rrig = <^5 = 0. In the presence of the effective /r-term we can neglect the 
original /i parameter, eliminating yet additional independent scale in the Lagrangian. Correspondingly, /i will 
in the following indicate the effective /i term. The relevant NMSSM Lagrangian now reads 

- 4 NMSSM ^ ^2 | 5|2 ^ (^xAxH^HdS + I + h.c.) . (3) 


_/-NMSSM _ _2 


“'soft 


= m^s l^r + + 


m 


'2 


In the MSSM, the minimization conditions of the Higgs potential can be used to replace and m|^^ by mz 
and tan (3 in the broken phase. Using the additional minimization condition of the NMSSM m| can be expressed 
in terms of /i. The Higgs-singlet sector m is therefore fully described by the parameters A, /^, Aa, A^, tan 
and the mass of the Z-boson. 

Eor specific NMSSM models we have to define the input scale of these parameters. The ratio of the VEVs 
tau/d is always evaluated at the weak scale m^, because it assumes electroweak symmetry breaking. Eor the 
high-scale models discussed in Sec. [ITT] A,/t: and /i are set at the SUSY scale of 1 TeV, while A^ and A^ can 
either be unified to Aq at the GUT scale or set individually (also at the GUT scale). Eor the low-scale models 


in Sec. [TT] all supersymmetric parameters including A, 
at the SUSY scale. 


Aa, A^,/i, the squark and slepton masses, etc. are set 


B. Higgs—singlet—singlino sector 

Compared to the minimal supersymmetric Higgs sector of the MSSM, the phenomenology of the NMSSM is 
strongly modified by the additional particles, a scalar and a pseudoscalar Higgs bosons and a fifth neutralino. 
While in general the mass of the singlet states is a free parameter, we will assume that the singlino contributes 
to a light LSP and that the singlet Higgs states are therefore lighter than their SM-like counterparts. Eor 
example, the SM-like Higgs boson with its mass of 125 GeV will typically be the second-lightest CP-even Higgs 
scalar. Eor a scale-invariant superpotential we can write out the symmetric Higgs mass matrix [33[ |34] in the 
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(i^, 5')-basis, where h is the SM-like Higgs boson 


^'^2/3 - 2 1 ^-^ 

^ V r y <52/3^1 


M; 


H,h,S 


= run 


(a 

C2/352/3 ^ 

2 2 
^2/3 +<^2/3^ 


■^2/3 (^a + 1 ^ 1 ^) 


2 A 


gmz 

/i-52/3 — 


S2/3/^/i 


gmz 

mz 


2g‘^li 


( 4 ) 


using the usual notation 5/3 = sin^^ and = cos p. In terms of the ordered mass eigenstates H 2 = i^i 25 means 
that throughout our analysis we identify the second-lightest Higgs with the observed SM-like state. Instead of 
n and A, the modified parameter set 

"=A 
A 

9 

appears in the diagonal entries for the light and heavy MSSM-like Higgs states. In the following, we replace n 
with R but keep A instead of trivially rescaling it by a constant g. At tree level the two NMSSM parameters 
take the pressure off the stop sector for small values of tan /3. In our basis conventions the second Higgs state 
it the SM-like observed resonance. This means we can decouple the singlet contributions from the observed 
Higgs. Setting 


(singlet mass parameter) 

(singlet decoupling parameter) (5) 


(4-A ( 6 ) 

Vs2/3 J 

removes the (2,3) entry from the mass matrix and therefore decouples the singlet sector from the SM-like Higgs 
boson . Note that this condition does not require any of the couplings in the NMSSM potential of Eq.(|^ to 
vanish. 

More generally, we can decouple the singlet from all other Higgs states in the limit A ^ ^ < 1. This way the 
two corresponding entries in the extended Higgs mass matrix vanish. To make the singlet itself heavy we need 
to increase its entry in the mass matrix in the limit X g. Neglecting the singlet entry in the Higgs mass 
matrix is ( 2 k/i)^, which for finite k consistently decouples with the single condition A ^ 1 . 


Aside from the extra CP-even Higgs, the singlet extensions of the MSSM Higgs sector adds an additional 
pseudoscalar. We can transform the weak eigen-basis {Hu, Hd, S) into a mass basis (A, S) by a rotation, so that 
A = cpHu + spHd. For large values of tan/3 the mass eigenstate A is approximately given by Hd- Removing 
the Goldstone modes the 3x3 mass matrix in terms of the weak eigenstates can be reduced to a 2 x 2 mass 
matrix in the basis (A, S'), which reads 


^A,S — ^ 


/ 2 /i (Aa + Rg) 
S2gm% 


gmz 

A^ /A 


(Aa - 2Rg) 


( 4/i2 


^2/3yy ( ^+2^) -3k 






rrin 


X n \ 


sl^ml 


gmz S2I3 

g "^z 


(7) 


ml / 


The pseudoscalar mass eigenstates are denoted as Ai and A 2 . In the second form we use the singlet decoupling 
condition Eq. (H and assume 52/3 ^1. As for the scalar sector, the limit X g decouples the singlet; its 
squared mass is then given by —SRgAu. The upper left entry of the matrix then corresponds to the MSSM 
pseudoscalar mass, m\ = 2/i (Aa + Rg) /<S 2 / 3 . This way we can choose either this MSSM-like mass or Aa as 
input parameter. Similarly, we can replace A^ with the lower-right entry in M\ ^ as the input parameter. 
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Finally, the supersymmetric partner of the singlet field, the singlino appears in the neutralino mass matrix. 


Ml 

0 

^^Z C/3 


0 

0 

M 2 



0 



0 


A 

C/3 


-9 

-mzsp- 




0 

A 

-mzc^- 

9 

0 

0 

A 

A 

-mzS/3- 

-mzC/3- 

2 k/i 



9 

9 



The bottom-right entry indicates that in accordance with Eq.([^ the combination 2^/i determines the singlino 
mass. The gauginos do not mix with the singlino. To altogether decouple the singlino, we have to remove the 
singlino-Higgsino mixing via A ^ 1 and at the same time make the singlino heavy, k ^ 1. In contrast, for 
^ < 1/2 the LSP will be mostly singlino. In this case the LSP mass is approximately given by the solution 
to 




2 k/i = 


rrin 


- I1S2I3 

g2 ^2 _ g2 


m% 

rUv — 2/^/i H—^ ^ 

^ 9^ ^ 


2 k - S2P 

4k^ — 1 


so that the LSP mass can be fixed by adjusting k. 


(9) 


For the interpretation of the galactic center excess a light pseudoscalar will be crucial. To describe its relevant 
couplings we have to rely on the different mixing matrices. The neutralino mass matrix will be rotated into 
its mass eigenstates through a matrix (Nij) with i^j = 1...5. To rotate the pseudoscalar mass matrix into its 
mass eigenstates we also have to consider the Goldstone mode. The corresponding mixing matrix is (Pij) with 
i = 1 , 2 and j = 1 , 2 , 3 because the Goldstone is not counted as part of the mass eigenstates Ai^ 2 - The lighter 
pseudoscalar Yukawa coupling to bottom quarks is given by 

_ • ^6 ^ /inA 

QAibb — Til 5 (^0) 

where Pn is the component of the lightest mass eigenstate. The coupling to the down-type squarks and 
sleptons is enhanced by 1/ cos /d, which is relevant for large tan (3. The coupling mediating the light pseudoscalar 
decay into the lightest neutralinos is given by 


QAixx —(Pii Y14Y15 + Pi 2 ^ 13 T /^15 + P13XI13N14,) — XkV^PisNi^ 

— — 92N12) (PiiY i3 — P12Y14) . (11) 

When we assign the second Higgs to be SM-like the lightest pseudoscalar will be mainly singlet. Therefore the 
coupling simplifies to 


9 Aixx — {XI13N14 — kNi^^ , ( 12 ) 

where we set Pn = P 12 ^ P 13 — 1. As Y 13 and Y 14 differ in sign, both contributions will add up. For a 
singlino LSP we have sizeable Y 15 ^ 1 , but following Eq.(|^ k ranges around m^/( 2 ya). This means that the 
singlino term in 9Aixx decreases with increasing but the same is true for the Higgsino fractions ^"13 and ^" 14 . 
Altogether, a large mediator coupling to the LSP points to a singlino LSP. 


C. Dark matter annihilation 

The lightest neutralino being the lightest supersymmetric particle of the NMSSM provides an excellent dark 
matter candidate. The key constraint for our parameter study will be its relic density assuming thermal 
production, leaving ns with distinct choices of annihilation mechanisms [lain]. Eirst, with a dark matter mass 
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around 40 GeV typical co-annihilation channels [TSHU with particles coupling to the Z-boson are constrained 
by Z-pole measurements at LEP. More generally, they are excluded by direct LEP searches, unless their masses 
are carefully tuned to only produce soft particles in their production and decays. What is left is neutralino- 
neutralino annihilation through an s-channel mediator: 


• Z-funnel annihilation through a Higgsino component in the interaction 

^ 9 ft [-^13-^13 - -^14^^14] • (13) 

Z cos uw 

Because of the velocity dependence of the annihilation rate {av) this channel usually prefers ESP masses 
slightly above or below 45 GeV, because directly on the Z-pole the annihilation is too efficient. 


• scalar iEi 25 -funnel annihilation, where the mass of the LSP has to strictly speaking be around 63 GeV. 
However, in combination with other annihilation channels the iEi 25 -funnel can give the largest contribution 
already for LSP masses around 55 GeV. The coupling to the Higgs can be found in Eq.(22). The Hxx 
coupling will be relevant for invisible Higgs decays, which is why we will discuss it in detail in Sec. im 

• pseudoscalar H-funnel annihilation. Unlike in the MSSM we now have a singlet and a Higgs pseudoscalar 
channel. Both of them can lead to a highly efficient annihilation with av (x . If the pseudoscalar is 
mainly singlet, the relevant contributions to the neutralino coupling in Eq.(ll) reduces to 


QAxx = (-^ 13^^14 - 


(14) 


The Higgsino components differ in sign, so that the absolute value of the couplings add up. Large Higgsino 
and singlino components therefore lead to a strong coupling. 


D. Data and tools 

In addition to the dark matter relic density and the observed Higgs mass value there exists a wealth of 
measurements which might constrain supersymmetric models. A large part of the data which we use for our 
SEitter parameter study are listed in Tab. [Tj 

The invisible width for Z-decays, P^^inv, is identified with the additional contribution to the Z-width from 
decays into a pair of LSPs. Eor LSP candidates with a mass smaller than 45 GeV, the LEP results for the Z- 
width can be powerful constraints. In addition, we require the lightest chargino to have a mass above 103 GeV, 
because it is very hard to avoid the LEP2 constraints for charged particle production [22]. This constraint 
becomes important when we consider regions with small /i < 200 GeV. If the sum of the lightest GP-even and 
CP-odd Higgs is smaller than the mass of the Z-boson, the total width gets an additional contribution. This 
contribution is compared to the difference between the SM prediction and the measured total width of the 


measurement 

value and error 


TflHi25 


(125.09 ± 0.21stat ± O.llsyst ± 3.0theo) GcV 

issrw 



0.1188 ± O.OOlOstat ± 0.0120theo 

m 



(287 ± 63exp ± 49sm ± 20theo) • 10“^^ 

m 

BR(B 

Xsl) 

(3.43±0.21stat±0.07syst)-10-‘‘ 

EH] 

BR(B° - 


) (3.2 ± 1.48tat ± O.Ssyst ± 0.2theo) ’ 10“® 

m 

BR(B+ - 

9 T+v) 

(1.41±0.43stat)-10-^ 

m 

Am^o 


(0.510 ± 0.004stat ± 0.003sy8t ± 0.400theo) • 10^2 
(17.69 ± O.OSstat ± 7.00theo) • 10^2 

m 

Am^o 


m 

P Z—>-inv 


(-1.9±1.58tat±0.2theo) MeV 

m 

Pz^Higgs 


(6.5 ± 2.3stat ± l.Otheo) MeV 

m 

mt 


(173.5 ± 0.68tat ± O.Ssyst) GeV 

m 

m + 

Xi 


> 103 GeV 

m 

(JxN 


direct detection limits 

m 


Table I: Data used for the fit including their systematic and statistical errors from the measurements and theoretical 
errors for SUSY calculations as far as they are considered. 
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Z. The mass of the top is an input parameter to the supersymmetric SM-like Higgs. The measurement of 
can only be satisfied with small slepton masses around 400 GeV or lower. As we decouple the sfermion sector 
this measurement will only lead to an overall constant contribution to the log-likelihood. Similarly, with the 
decoupled stop sector the exact value of the SM-like Higgs mass is irrelevant to our analysis, because it can 
always be re-adjusted independently of our parameters of interest [34]. Finally we also include the XenonlOO [43] 
limits on direct detection. They are most powerful in the 30 to 100 GeV range, so this measurement will prove 
to have a strong exclusion power for the scenarios we are interested in. 

The Hooperon as an explanation of the galactic center excess |4| is described by two parameters: the LSP mass 
and the annihilation cross section in the center of the galaxy. For our analysis we will assume the conservative 
range m 


=(30 ■ ■ ■ 70) GeV 


av =(0.4---10) ■ 10“^® ^ 


(15) 


A more detailed analysis in the two-dimensional plane will be part of the full analysis in Sec. HB 


The SFitter fit then determines multidimensional likelihood maps for the model parameter space. A set of 
Markov chains selects points in the model space following a Breit-Wigner proposal function. For each point we 
compute all considered observables and determine a generalized value [441446] . Theoretical and experimental 
errors are combined using the RFit scheme im, which corresponds to a consistent profiling over the nuisance 
parameters. Gorrelations between observables like different channels at the LHG are taken into account. These 
likelihood maps can be projected on two-dimensional planes as well as single-parameter measurements using 
profile likelihood or Bayesian methods. Throughout this paper we will apply frequentist profile likelihoods and 
this way avoid volume effects. 

For the TeV-scale fits the same technique is used to select points of the Markov chains following the proposal 
function. When we display observables instead of the likelihood function, only points that pass sharp criteria, 
e.g. Vtyji? within the theory uncertainties, are displayed. Similar to a profile likelihood, we assign the value of 
the point with highest likelihood to each bin. 


For constraints specifically for the NMSSM we rely on NMSSMTools [48] to calculate the supersymmetric 
mass spectrum, the Higgs branching ratios, the B observables, {g — 2)^, and electroweak precision observables. 
The relic density and the direct detection limits are calculated with MicrOMEGAs [49]. The number of 
events in the Higgs production channels at the LHC for the SM-like i7i25 are computed using the standard 
SFiTTER-Higgs setup [50], HDecay [51], and NMSSMTools. 

For the MSSM study in the appendix the supersymmetric spectrum is calculated using SuSpect3 [52] while 
the Higgs branching ratios are computed using Susy- Hit[5T] and HDecay. For the electroweak precision 
observables we rely on SusyPope [53]. Finally, we use SuSpect3 to compute the H-observables and {g — 2)^. 
As for the NMSSM the relic density and the direct detection limits are calculated with MicrOMEGAs. 


II. TEV-SCALE NMSSM 

The NMSSM can explain or motivate interesting measurements which are not accessible in the MSSM. This 
is typically an effect of the extended Higgs particle content, for example including a second pseudoscalar state. 
One of these measurements is the galactic center excess in gamma rays, based on light neutralinos and a light 
pseudoscalar mediator. Another is invisible decays of the SM-like Higgs boson. It will turn out that both of 
them are strongly linked. 

We assume a Zs-symmetric NMSSM, where all input parameters are set at the scale 1 TeV. The absence 
of unifying assumptions leads to a large number of model parameters, namely the slepton and squark masses, 
the trilinear couplings and the gaugino mass parameters. In this study we decouple the squarks, sleptons, and 
gluinos by setting the soft masses to 10 TeV and the trilinear couplings to zero, because these sectors are not 
experimentally relevant. To obtain the correct 125 GeV Higgs mass we adjust the stop masses in the TeV-range 
appropriately. For the Higgs-singlet sector the relevant input parameters are 




( 16 ) 
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Alternatively, we can replace Ax^k by the diagonal entries in the pseudoscalar mass matrix of Eq. 0- The 
neutralino-chargino adds the free parameters Mi and M 2 . 


While is given at the SUSY scale of 1 TeV, NMSSMTools computes the Higgs masses at the averaged 
squark masses, in our case around 10 TeV. In Sec. [ITT] we will discuss the scale dependence of the singlet related 
parameters in detail. The approximate form of A^ at the averaged squark mass scale can be computed from 
the RGE via 


A^{10 TeV) = A«,olO^ + f 10^ - ij 


4 31ogl0,^ 

^-^,0 + 4^2 • 


(17) 


Due to the large value of |Aa|, that further increases with the absolute value of /r, A^ increases from —250 GeV 
to the order of 30 GeV at 10 TeV for A = 0.2 and Ax = 8 TeV. When we consider small singlet masses, the 
scale dependence of A^ plays an important role in their parameter dependence. 


A. Galactic center excess 


With the galactic center excess of gamma ray photons we add an experimental motivation for a light neutralino 
in combination with a light pseudoscalar mediator mn] to the (largely negative) experimental results in Tab|^ 
The spherically symmetric excess can be explained by annihilating dark matter, more specifically by a neutralino 
in the range of 30 GeV to 40 GeV or even 70 GeV, that annihilates into a pair of fermions, for example XX 
Eor a type-II two-Higgs doublet model the bottom and tan Yukawa couplings are aligned, which means that 
the assumed decay to a bb pair will dominate over the accompanying decay to . In our preferred data 

regions the Ai decays into bb (94%) as well as into rr (6%) pairs, the ratio reflecting the size of the Yukawa 
couplings and the color factor in the case of quarks. 

In the MSSM such a light neutralino ESP is a problem, because it will lead to too large a relic density. In 
the presence of direct SUSY and Higgs search results there is no obvious way to annihilate it efficiently enough 
to arrive at the observed relic density. Different co-annihilation channels |T9]-[2T] require new charged states in 
reach of LEP2 and are therefore not viable [22] . 

In contrast, the two pseudoscalars in the NMSSM can mediate a sufficiently fast annihilation, because the 
LSP annihilation through its resonance pole is not p-wave suppressed like it is for vector bosons or scalars. 
While the mass of the MSSM-like pseudoscalar Higgs is strongly constrained by heavy neutral and charged 
Higgs searches, the additional light pseudoscalar can be mostly singlet. The neutralino coupling to such a 
mediator is given in Eq. 0 for a largely singlino mediator turning into Eq.([l4|). 


We consider a generic NMSSM scenario based on a light singlino with a bino and a Higgsino admixture, 
f.e. k ^ 1 [27]. An alternative solution presented in Ref. [30] combines a bino-Higgsino LSP with an NMSSM 
pseudoscalar mediator, but will not lead to the measurable invisible Higgs decay we are interested in. Light 
neutralinos with a sizeable wino or Higgsino component are essentially ruled out by Z-pole measurements and 
by chargino searches at LEP, so we decouple the wino at M 2 = 3 TeV. In this limit we can link k to the 
neutralino mass through Eq. 0. We then adjust Ax such that the singlet component of the SM-like Higgs is 
minimized, avoiding Higgs sector constraints altogether m- 

If the annihilation process leading to the galactic center excess proceeds via a pseudoscalar decaying into a 
pair of bottom quarks, today’s dark matter annihilation cross section is [24] 


(TV 


v=0 



^Mxx 


9A,bb'rn\ 


X (my-4m|) +m\Tl 


(18) 




Large tan^d increases the coupling to the bottom quarks given in Eq.(lO). A strong coupling to the bottom 
quarks also leads to a small width of the light pseudoscalar which further increases the cross section when one 
is close to the on-shell condition. 
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X 

Figure 1: Neutralino relic density (left) and mass of the light pseudoscalar Ai (right) color-coded as a function of // 
and A. The two lower panels are zoomed into the respective upper panels. The orange regions in the lower left panel 
are compatible with the relic density = 0.107 ... 0.131, considering the theoretical uncertainty. In addition to 

the usual decoupling through large scalar masses we hx tan/3 = 40, = —250 GeV, according to the decoupling 

condition Eq. and k corresponding to an LSP mass of 40 GeV through Eq. §. 


Using Sfitter we analyze the NMSSM parameter space with a focus on the SM-like Higgs mass, the LSP 
mass, the observed relic density the Z width measurements, the Xenon direct detection constraints, 

and the galactic center excess. As expected from Sec. |IB| the key model parameters are the Higgsino mass 
parameter /i, the singlet mass parameter k, and the coupling A, which links the singlet to the MSSM Higgs 
sector. In Fig. we observe a broad band in the A vs /i plane, which is defined by non-zero values for the Higgs 
masses: for large /i ^ ruz and A ^ 1 the mass of the lightest CP-odd scalar Ai vanishes. This follows from 
Eq.(pT|) and the leading term in Eq.(|^, such that the pseudoscalar mass matrix in Eq. 0 only depends on /i 
and A. A^ is of the order of 100 GeV at the scale of the averaged squark masses. Following Eq. 0 , an increase 
of A and Ay, leads to an increase of A^. This results in a smaller mass eigenvalue of Ai. As Ax increases with 
/i following Eq. 0, higher values of /i lead to smaller pseudoscalar mass. 

The mass of the lightest CP-even scalar Hi is given by the (3, 3) entry of the Higgs mass matrix in Eq.Q in 
the Higgs decoupling limit Eq.([^, 


^H,h,S 


= rrir. 


33 


A^ 


^ (1 — 52/3^) H-^ [A,^ + 4k/i) 




(19) 


enters here with a positive sign, so that for small /i < mz and A ^ 0 the mass of this CP-even singlet 
vanishes. Possible experimental constraints are expected to further reduce this band. 


Within this broad band shown in the upper panels of Fig. [^the structure originates from two sets of input 
parameters to the calculation of the relic density. On the one hand there is a strong dependence on the mass 
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of the pseudoscalar mediator, on the other hand the couplings of the LSP depend on the gaugino and Higgsino 
content of the LSP. 

In the upper part of the band, with rriAi < 80 GeV and 0.26 < A < 0.3 the dark matter annihilation is 
mediated by the Z-funnel, with a coupling to the Higgsino content proportional to as given in 

Eq.(13). Smaller values of A, correlated with large values of /i, decouple the singlet/singlino sector from the 
MSSM. An efficient dark matter annihilation is not possible, and the relic density is too large. On the other 
hand, too small values of /i and large A increase the Higgsino-singlino mixing via the off-diagonal terms in the 
neutralino mass matrix. The Higgsino component in the relic neutralino then results in too small a relic density. 
In between. Fig. indicates a well-defined regime with the correct relic density. The corresponding mass of the 
lightest pseudoscalar Ai indicates that the resonance condition ~ 2m^ is only fulfilled at the lower end of 
this regime, while the larger part of the allowed band relies on Z-mediated dark matter annihilation. 

In the lower part of the band, defined by the onset of the resonance condition ~ a steep decrease 
of the relic density leaves a very narrow strip where the annihilation proceeds via the A-funnel and reproduces 
the correct value of With the increasing A-funnel contribution, the Z-mediated annihilation rate has to 

decrease, which means that the allowed region bends towards larger values of ji. The moment the resonance 
condition is actually fulfilled, the annihilation through the A-funnel becomes too efficient, and the predicted relic 
density drops well below the measured value. For tyia^ ~ 90 GeV, corresponding A < 0.24 and /i ^ 300 GeV, 
the annihilation again proceeds off-shell, predicting the correct relic density starting with a reduced Z-mediated 
rate at large ja. At A ~ 0.225 and /i ~ 275 GeV the annihilation is again mediated by the Z-boson alone. 

A few hardly visible points with the correct relic density at the very top of the allowed mass band arise from 
iLi-funnel annihilation and will be of no relevance to our further discussion, because the scalar mediator with 
its p-wave suppression fails to explain the galactic center excess. 


The lower panels of Fig. focus on Ai-funnel annihilation just below the resonance condition. There are 
two regions where the annihilation cross section is compatible with the galactic center excess — within a 
comparably broad range of av = (0.4 ... 10) • 10“^^cm^/s. The region below the resonance condition rriAi ^ 
is compatible with the relic density, while the other one is not. 


B. Invisible Higgs decays 


Invisible Higgs decays have long been in the focus of LHC studies BE]- At the LHC the upper limits 
on an invisible branching ratio are 57% in the WBF channel from CMS [8] and 78% combining associated 
Higgs production and gluon fusion from ATLAS [7]. Eventually, the HL-LHC should be sensitive to invisible 
branching ratios of a few per-cent in the WBF production channel [2]. Usually, such invisible Higgs decays 
are associated for example with a Higgs portal to a scalar dark matter sector [32]. We will show that in the 
NMSSM, invisible Higgs searches can also probe a relevant part of the dark-matter-related parameter space 
through the decay i7i25 ^ XX- Because this decay requires relatively light LSP neutralinos these scenarios can 
be linked to the galactic center excess discussed in Sec. H A[ 


The decay width of the CP-even Higgs into two neutralinos is given by 


r(ifi25 ^ xx) = 


with A = 1 — 


4m| 


m 


( 20 ) 


H 


The Higgs-LSP coupling in the MSSM is 


9hxx 


= (^^11 — g'N12) (sinaVis + cosaVi4) 


MSSM 

= (pVn - g'N^2) (^2iA^i 3 + ^22A^i4) • 


( 21 ) 


The Nij are the entries of the neutralino mixing matrix, and S 2 i are the entries of the CP-even Higgs mixing 
matrix. In the simple (2 x 2) case the latter are expressed in terms of the mixing angle a. In the MSSM 
invisible Higgs decays have to be mediated by gaugino and Higgsino fractions combined, or more specifically 









11 




|i=190 GeV 



Figure 2: LSP mass, mass splitting between the LSP and the pseudoscalar mediator, and invisible Higgs branching ratio 
with /i = 190 GeV and = —250 GeV in the k — Mi, k — X, and k — Mi planes. All displayed points are compatible 
with the relic density. Xenon, a chargino mass above 103 GeV, and the correct SM-like Higgs mass. Moreover, they 
always have an invisible branching ratio BR(Mi 25 ^ xx) > 10%. The black points are consistent with the galactic 
center excess. 


by a mixed bino-Higgsino LSP. In the NMSSM this coupling receives additional contributions from the singlet, 
namely 


9hxx = 9hxx 


+ ^/ 2 A [(^21 Vi4 + 522Vi3)iVi5 + 523iVi37Vi4 - ^ 523 ^ 25 ] 


MSSM 


( 22 ) 


Now, invisible Higgs decays can proceed to bino-Higgsino, singlino-Higgsino, and in the presence of a singlet 
component in the Higgs boson even pure Higgsino and pure singlino LSPs. 


Just like the galactic center excess, invisible decays of the SM-like Higgs benefit from a light, mixed neutralino 
LSP. Decoupling squarks, sleptons, and gluino we can ask if there are regions with a large branching ratio 
^125 ^ XX foi* a mixed bino-Higgsino-singlino LSP, and such parameter regions can be related to the galactic 
center excess. 

Following Eq. § the mass of the singlino LSP is given by 2/i^. To suppress the Higgsino component we 
require ^ < 1/2 and keep /i > 190 GeV in a first step. This way the Higgsino component in the LSP ranges 
around 5% to 10%, leading to a sizeable coupling to the Higgs sector but preventing a large coupling to the 
Z-boson. Our choice of parameters for the galactic center excess fixes tan/3 = 40 and < —250 GeV. 
The coupling condition Eq. (H determines Ax to minimize the singlet component of the SM-like Higgs. The 
remaining parameters are A, k and Mi. Small values of Mi increase the bino component of the LSP while small 
values of k increase the singlino component. To keep a constant hooperon mass, k and Mi have to behave 
inversely proportional. 


In Eig. [^we show the result of the SEitter analysis, starting with fixed fi = 190 GeV. Of the experimental 
constraints discussed in Sec. |ID| we now only include the relic density, the Xenon constraints, the chargino 
mass constraints, the invisible Z-width constraint F^^inv < 2 MeV and the SM-like Higgs mass constraint 
rriH = 122 ... 128 GeV. We only show parameter points with BR(iLi 25 ^ xx) > 10%. Eor our starting value 
/i = 190 GeV we fix A,^ = —250 GeV. 

Three distinct, narrow strips for example in the upper row of Eig. are defined by constant LSP masses 
around 40 GeV, 48 GeV, and 55 GeV. Eor ~ 55 GeV the annihilation is mediated by H 125 . If the mass 
moves closer to the Hi 2 s on-shell condition the relic density becomes too small. At the lower end of the 55 GeV 
strip the additional annihilation through the Ai pseudoscalar and the Z becomes too strong to reproduce the 
observed relic density. Numerically, the reason is that oc k reaches 120 GeV around k = 0.16, which opens 
a pseudoscalar-mediated LSP annihilation channel. Eor /i = 190 GeV this coincides with the possibility to 
efficiently annihilate though an s-channel Z-exchange. 
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The 40 GeV and 48 GeV strips are defined by the Z-mediated LSP annihilation. Each of them lives on one 
side of the Z-pole, because the annihilation on the pole is too efficient to give the correct relic density. Both 
strips follow the asymptotic behavior of the LSP mass. The 40 GeV strip continues towards high values of 
Ml, but with a reduced LSP mass. The reason is that the additional annihilation mediated by Ai adds to 
Z-mediated annihilation. Finally, the annihilation via the pseudoscalar connects the annihilation channels for 
= 48 GeV and = 55 GeV around Mi = 70. 

In the broad regions with Mi = 10 ... 40 GeV the annihilation is in addition mediated by Hi {k > 0.1) and 
by a combination of Hi and Ai (^ < 0 . 1 ). 

One key feature is the hole in the allowed parameter space around k = 0.15 and Mi = 50 GeV. For example 
along the 40 GeV strip the LSP composition changes from bino to singlino, with a 5% to 12 % Higgsino 
contribution. While the sum of the two Higgsino components increases towards the singlino LSP, their ratio 
switches. This leads to an intermediate region where both components have a similar value. At this point gzxx 
given by Eq.(13) vanishes, interrupting the Z-mediated annihilation. The other annihilation channels are weak, 
so the relic density is too large. 


In the right panel of Fig. we see that the parameter points which can account for the galactic center 
excess appear precisely where we also expect invisible decays for the SM-like Higgs H 125 . Indeed, for Mi = 
50 ... 60 GeV and k ~ 0.1 we find the highest branching ratio Hi 2 s XX- This region has unique properties: 
as discussed above, the dark matter annihilation proceeds through an on-shell Z-boson, with = 45 ± 2 GeV. 
The LSP is a mixed state with 8 % Higgsino, 30% to 50% singlino, and 40% to 50% bino content. The two 
Higgsino components are of the same size, strongly reducing the Z-coupling gzxx' lightest pseudoscalar Ai 
has a mass around 115 to 135 GeV, but as an at least 95% singlet remains undetected. Because of the strongly 
mixed LSP content, each Nn term can contribute to gH 2 XX' large value induces an invisible branching 

ratio of up to*. 


BR(Mi 25 ^ XX) < 70% for g = 190 GeV , (23) 

for the parameter space consistent with the galactic center excess. 

In this region with an overwhelming invisible Higgs decay width we need check a few experimental constraints 
not represented in Fig.|^ First, the cancellation of the two Higgsino components renders the partial decay width 
Pz^inv smaller than 0.1 MeV. For a mostly-singlino LSP region it increases through large V 14 values up to 
2 MeV. The Xenon limits on direct dark matter detection exclude the region centered around Mi = 40 GeV 
and k = 0.2, linked to the Hi and Ai annihilation channels. However, none of these additional constraints 
affect the parameter regions with invisible Higgs branching ratios between 10% and 30%. 

In the next step we vary /i, and with it A^ to keep the singlet Higgs masses in the desired range. This means 
that for g = 220 GeV and g = 320 GeV we have to set A^ = —280 GeV and A^ = —350 GeV, respectively. 

Increasing g slowly impacts the LSP annihilation channels. First, we see that the allowed regions for g = 
220 GeV shown in the upper left panel of Fig.j^are very similar the case of g = 190 GeV. This indicates that our 
above results are not very fine-tuned. For a fixed LSP mass an increase of g merely leads to a smaller Higgsino 
component, which in turn leads to a smaller Zxx coupling. It also decreases the invisible Higgs branching to 

BR(Mi 25 ^ XX) < 40% for g = 220 GeV , (24) 

in the relevant parameter region for the galactic center excess. 

For larger g = 320 GeV the picture changes: for a fixed LSP mass an increase of g leads to a smaller Higgsino 
content. While for g = 190 GeV the two Higgsino components add to 5% ... 17%, they now stay in the 2 % ... 5% 
range. This immediately leads to a smaller Zxx coupling — and a reduced invisible Higgs branching ratio. 
The Zxx coupling implies that the LSP mass has to be closer to the on-shell condition to give the correct relic 
density. Indeed, for a bino-like LSP and g = 320 GeV, the lowest mass strip is now defined by ~ 42 GeV 


In this discussion we ignore the limits on invisible Higgs decays for example from the SFitter Higgs analysis There, we 

find BRinv < 30.6% at 95% C.L. in an 8-parameter coupling analysis. For a dedicated NMSSM fit this limit is expected to be 
even more constraining. 
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|i=220 GeV 



11=320 GeV 



Figure 3: LSP mass, mass splitting between the LSP and the pseudoscalar mediator, and invisible Higgs branching ratio 
now for 11 = 220 GeV (upper row) and fi = 320 GeV (lower row) in the k — Mi, k — X, and k — Mi plane. Gorrespondingly 
we choose = —280 GeV and -380 GeV. As in Fig.j^all displayed points are compatible with the relic density. Xenon, 
a chargino mass above 103 GeV, and the correct SM-like Higgs mass. Moreover, they always have an invisible branching 
ratio BR(iVi 25 —> xx) > 10%. The black points are consistent with the galactic center excess. 


instead of 40 GeV. Similarly, the high-mass strip moves down from = 48 GeV to 46 GeV. Altogether, in 
the lower left panel of Fig. we see that the annihilation regions mediated by the Z-fnnnel and the Hi funnel 
with = 55 GeV clearly separate. 

The annihilation processes can now best be identified in the central lower panel of Fig. showing the 
correlation between k and A. Annihilation through a Z-boson occurs in the two parallel strips with k ~ 0.065 
and k ~ 0.07. They are divided by the actual on-shell condition 2m^ = for which LSP annihilation 
becomes too efficient. 

Following Eq.Q and replacing Ax through Eq.(|^ we see that the Hi mass increases with A directly, as well 
as indirectly via A^. Larger values of A lead to a steeper increase of the scale dependence and thereby increases 
A^ at the 10 TeV scale. At the same time the neutralino mass increases with k. This explains why for an 
annihilation via Hi we find a strip increasing from k ~ 0.02 and A ~ 0.24 to ^ ~ 0.06 and A = 0.4 

Eor the galactic center excess the annihilation via the pseudoscalar Ai is crucial. The LSP mass decreases 
with K, while the mediator mass rriAi decreases with A. This is caused by the renormalization group running, 
which for larger A pushes A^ to larger values at 10 TeV. Eollowing Eq. 0 the pseudoscalar mass includes a 
factor —A^^ which means it indeed decreases with increasing A. To maintain the relation between the LSP and 
mediator masses, A and k have to be anti-correlated. This is what we observe in the two Ai-mediated strips in 
the central lower panel of Eig. These strips with the efficient pseudoscalar mediator also accommodate the 
galactic center excess, as expected from the simplified model analysis [24]. Again, as before the Ai-mediated 
annihilation blends in with Z-mediated annihilation. 

Eor 11 = 320 GeV we find an increasingly small number of parameter points which accommodate the galactic 
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Figure 4: Correlation between the LSP mass and the dark matter annihilation rate for our three reference points in /x, 
to be compared to for example Fig. 3 in Ref. m- The lines correspond to the two sigma limits there. If available, we 
indicate the leading LSP annihilation channels. Points with < 30 GeV shown in this figure are not included in our 
analysis and in Figs. and 


center excess as well as the current relic density. However, they reside in a regime with 

BR(iLi 25 ^ XX) < 15% for /x = 320 GeV , (25) 

again consistent with the galactic center excess. One aspect which will become relevant when we embed the 
NMSSM in a unified framework at the GUT scale is the Mi-dependence of the Hooperon-compatible parameter 
points. For all three /x-values we observe a tail towards large Mi values, for which the LSP properties do not 
change any longer. While we do not show values above Mi > 90 GeV, we could extend the curves to much 
larger bino masses. 

For a combination of the different /x values assumed in the above analysis we show the two-dimensional corre¬ 
lations between the LSP mass and the annihilation rate in Fig.|^ The valid NMSSM parameter points include 
LSP masses below 30 GeV, which we do not consider in our analysis of the invisible branching fractions mn]. 
We see that most of the NMSSM points would prefer a larger annihilation rate than the central Hooperon 
value, but this annihilation rate can be accommodated by moving the different masses slightly on and off the 
respective resonance conditions. Moreover, the majority of allowed points have LSP masses between 10 and 
30 GeV, in particular for /x = 220 GeV. On the other hand, for all three values of /x it is possible to enter the 
preferred region taken from Ref. [ 12 ] in the two-dimensional av vs plane. 


III. HIGH-SCALE NMSSM 

Instead of using the full set of TeV-scale model parameters we can require the NMSSM to fulfill a set of 
unification assumptions at large energy scales. We start with a unified squark and slepton mass mo, a unified 
gaugino mass mi/ 2 , and a unified trilinear coupling Aq at the GUT scale [55]. Furthermore, we require a Z 3 - 
symmetry to remove for example the /x-term and replace it with an effective /x-term induced by the singlet VEV. 
Because we do not require unified Higgs masses and m| we refer to the model as the NUH-NMSSM 

(non-universal Higgs masses) [56] . 

From Sec. m we know that to produce a light scalar and a light pseudoscalar, k has to be small and 
evaluated at 10 TeV has to range around the electroweak scale. The running of the different model parameters 
from the GUT scale to the weak scale is described by renormalization group equations, for example for the 
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singlet couplings 


IGtt^ 

IGtt^ 

IGtt^ 


dX 

d\ogQ‘^ 

dn 

d\ogQ‘^ 

dk 

d\ogQ‘^ 


— (4A^ + 3h^ + 3hl + 2 k,^ ~ 91 ~ ^qI) 

I (6A^+6,.^) 

— (4^^A^ + 2A^ — 3h^ — 3hl — ^ g\ ^ 3gl) . 


(26) 


The Yukawa couplings are defined as mj = hfV sin f3/^/2. The couplings k. and A appear squared, so that for a 
choice of signs our argument will hold for their absolute values. If we neglect the gauge couplings, A increases 
with Q and runs into a Landau pole. If we also neglect the Yukawas, which accelerate this increase, the Landau 
pole is approximately given by 


^{Q) — ^0 1 


Ag 1 Q 


- 1/2 


oo . 


(27) 


For our theory to be valid up to Q = 10^^ GeV the Higgs-singlet coupling is limited to Aq < 0.81 at Qo = 1 TeV. 
The large and also increasing top Yukawa coupling further accelerates the approach of a strongly interacting 
regime, requiring Aq < 0.5 ... 0.6 for a valid theory to the GUT scale. Assuming roughly constant A and also 
ignoring the Standard Model Yukawa and gauge couplings, the running singlet self-coupling k is given by 


f^{Q) = 



(28) 


The maximum value of k, for which the theory is defined up to Q = 10^^ GeV ranges around 


^max 




A =0 
A = 0.6 . 


(29) 


In this approximation we can compute a few example values by iterating: starting with A = 0.3 = /^o at the 
TeV scale we find k = 0.43 and A = 0.61 at 10^^ GeV. The singlet mass parameter k decreases from 1.0 to 0.7. 

In the running of k the top Yukawa coupling enters with a negative sign. Therefore k increases slower with 
the scale than A as long as /^,A ^ If we consider larger values as weak-scale starting points, k increases 
faster than A. A = 0.45 = kq gives k = 1.7 and A = 1.2 at the GUT scale, so k increases from 1.0 to 1.4. 


For the component fields both, the Higgs-singlet coupling and the singlet self-coupling come with associated 
mass scales. They run like 


dA\ 

IOtt^——= A}?Ax + 3h^At + 3hlAij + h^Aj- + 2k‘^A^ -h g\Mi + 3^|Af2 
d log 

167r^ (ik^g2 ^ + ^Hd + + 3k^m% -h k?A^^) 


(30) 


The increase of these mass scales towards high energy scales clearly does not help with the appearance of a 
strongly interacting Higgs-singlet sector in the NMSSM. 


Finally, we can ask how the new NMSSM Higgs-singlet parameters affect the running of the MSSM-like 
parameters. An interesting parameter in the MSSM is the stop mixing parameter, which now runs like 

IOtt^ ^log ~ + hlAi) + X^Ax + —glMi + 3^|M2 + —g^Ms (31) 
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While there will be an effect of the additional singlet on the running of the MSSM-like parameters, its impact 
will be numerically small. The only exception appears when we allow with Higgs-singlet sector to become 
strongly interacting at relatively low scales, in which case for example the stop mixing parameter will also 
sharply increase. 


A. Global analysis 


From the previous discussion it is clear that there are several, more or less distinct regions of the NMSSM 
parameter space, which allow us to describe the Hooperon. This is less clear when we constrain the model 
parameters through unification assumptions. Most notably, the unification of the gaugino masses at the GUT 
scale links the bino mass Mi to the gluino mass, which in turn is constrained by LHC searches 157]. 

Before we focus on the galactic center excess and invisible Higgs decays, it makes sense to test how well the 
unified NUH-NMSSM can accommodate all other available data listed in Tab. [Ij including the observed relic 
density and the Xenon limits on direct detection. We also include the SM-like Higgs couplings strengths from 
the SFiTTER-Higgs analysis [50j. Our parameters of interest are A, and ja. 


In the NUH-NMSSM with decoupled scalars (mo = 2 TeV) we face two major differences compared to the 
TeV-scale model. First, unification links the binos mass Mi to the gluino mass, which is constrained by direct 
LHC searches m- In addition, the SUSY-breaking singlet parameters Ax and are set at the GUT scale. 
This means that we cannot apply the simple decoupling relations for example given in Eq. (§. 

To get some control over the parameters we start with an SFitter likelihood analysis, for example to see 
how the requirement H 2 = M 125 translates into the high-scale parameters Ax, A^^mij 2 , and Aq. The unified 
gaugino mass is proportional to the gluino mass and hence constrained to be mi /2 ^ 500 GeV. The mass scales 
Ax and A^ will eventually be constrained by the requirement of a light scalar mass, and we fix their ranges 
to A^ = —1.5... 1.5 TeV and Ax = —1...5 TeV. From our experience with the TeV-scale NMSSM we limit the 
singlet parameters, still defined at 1 TeV, to A < 1 , /i: < 1 , and /i < 400 GeV. The ratio of the VEVs is fixed 
again to tan /3 = 40. 

As before, we use all data listed in Tab. |Tj In the Higgs sector we identify the observed SM-like Higgs 
with the second-lightest NMSSM Higgs and include the Higgs couplings measurements from ATLAS and CMS 
searches [50]. A set of two-dimensional profile likelihood projections is displayed in Eig. 


The upper left panel shows the profile likelihood projection on the k — X plane. We identify three distinct 
regions through their dark matter annihilation channels m- 

1 . a broad band with k = 0.1...0.3, A < 0.25, and fi = 90...200 GeV. It includes two LSP mass regions with 
an annihilation through Z- and iLi 25 -exchange. 

2 . a narrow strip around k ~ 0.42, A < 0 . 2 , and /i = 90...350 GeV. It relies on a light chargino either for 
co-annihilation or for t-channel exchange for efficient dark matter annihilation. 

3. a bulk region with k = 0.3...0.7, A > 0 . 2 , and ij. = 90...150 GeV. Here, the annihilation occurs through a 
mix of channels, notably including the light singlet pseudoscalar. 


The transition between the second and third region is not uniquely defined, but involves it the appearance of 
the Ai-funnel annihilation and a drop in the LSP singlino content from 70% ... 90% to 10% ... 70%. Hollowing 
Sec. HB[ an invisible branching ratio needs a LSP mass smaller than 62 GeV. Therefore the only region 
compatible with the galactic center excess and an invisible branching ratio will be k = 0.1 ... 0.25. 


The upper right panel in Eig. shows the profile likelihood projection on the k — plane. A lower bound 
/i > 90 GeV arises from the chargino mass limit. An upper bound is connected to the requirement that the 
second-lightest NMSSM Higgs be the SM-like state: the mass of the singlet-like Higgs is proportional to /i, so 
for large /i it approaches 125 GeV. This translates into the globally observed /i < 400 GeV. 

We then combine the range in /i with the mi /2 dependence. As mentioned before, the gluino bound sets a 
lower limit on mi /2 > 500 GeV. The combination of /i < 400 GeV and mi /2 > 500 GeV results in a sum of 
the bino and wino LSP components to be less than 1% throughout the plane. The mass and composition of 
the Higgsino-singlino LSP is set by /i, and A. Eor k > 0.5 it is mainly Higgsino, with its mass set by /i. Eor 
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Figure 5: Profile likelihood projections of the NUH-NMSSM assuming a SM-like H 2 = H 125 . Ah measurements shown 
in Tab. and the SM-like Higgs couplings m are included. Of the shown parameters and Ax are dehned at the 
GUT scale. The singlet parameters A, k and /x are dehned at 1 TeV, just like in the low scale scenario. 


k < 0.4 the LSP is mainly singlino, and following Eq.(|^ its mass is given by 2^/i . Large values of A lead to a 
stronger mixing between singlino and Higgsino. 

Of the list of regions introduced above we first consider the band with k = 0.1...0.3, where the singlino 
component is larger than 85%. The corresponding values of /x range from 90 GeV to 200 GeV, resulting in two 
regions of neutralino mass compatible with the measured relic density: for = 40 ... 50 GeV the annihilation 
is mediated by a Z-boson, while for = 55 ... 60 GeV the LSP annihilates via the SM-like Higgs into bb and 
partially into light Higgs bosons. Annihilation via the lightest pseudoscalar, as relevant for the galactic center 
excess can occur, but it is not a main annihilation channel. 

In the narrow second strip, the higher value of ^ = 0.42 leads to a smaller Higgsino component, that varies 
between 70% and 90%. The higher value of k leads to a mass ratio of about 0.85 between the LSP and the 
Higgsino-like chargino. This opens the chargino co-annihilation channel or neutralino annihilation into WW 
and ZZ through a chargino in the t-channel. In the upper right panel we can verify that this channel is open 
up to /X = 350 GeV. 

Finally, for A > 0.2 the different annihilation processes are no longer well separated. Inversely correlated to 
k = 0.3 ... 0.7 the singlino component decreases from 70% to 10%. For this region we find an upper bound 
of /X < 150 GeV, resulting in neutralino masses between 60 and 100 GeV. For neutralinos around 60 GeV the 
annihilation proceeds via the SM-like Higgs, while for larger masses the annihilation channel is a mixture of a 
pseudoscalar funnel, chargino co-annihilation, and t-channel annihilation via a chargino. 

The lower panels of Fig. [^show the profile likelihood projection onto the X — A^ and the Ax — planes. 
From Eq.Q and Eq.Q we know that the scalar and pseudoscalar singlet mass terms increase with A. The 
singlet-like scalar has to remain lighter than 125 GeV, leading to an upper limit on A depending on and 
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Ax- Large values of can either push the scalar mass to too large values or lead to a vanishing pseudoscalar 
mass. Both effects set an upper limit on A^. Because A^ is set at the GUT scale, small starting values of A^ 
can turn negative towards the TeV scale, leading to a very light scalar. The correlation between A and A^ 
occurs because for large A we need larger values of A^ to keep the pseudoscalar singlet heavy enough. 

In the right panel we see that large values of A^ are only possible for even larger values of ^4^. The scalar 
and the pseudoscalar singlet mass-squares differ by Am^ = 4:{A^kfi + neglecting the subleading term 

proportional to 171 ^ 2 :^ 213 ^^For both, the scalar and the pseudoscalar masses to be above zero, this mass 
difference cannot be larger than the actual mass scale. This means that large A^ has to be accompanied by 
even larger ^4^. 


B. Galactic center excess 


In the TeV-scale NMSSM a singlino-like LSP with a small Higgsino component can generate the galactic 
center excess in agreement with the relic density and linked to an enhanced branching ratio iLi 25 ^ XX- By 
definition, the NUH-NMSSM contains only a subset of the NMSSM models: the unification condition on mi /2 
impacts the range of Mi, and the stop masses can no longer be set independently of the remaining sfermion 
masses. 

To be consistent with the TeV-scale study we decouple the sfermion sector at mo = 10 TeV. This is compatible 
with the observed Higgs mass, when we adjust Ax accordingly. Again, tan/3 is set to 40, to provide a large 
coupling between the pseudoscalar Higgs and the down-type quarks. As before, /i is set at the SUSY scale of 
I TeV and limited to (150 ... 220) GeV which will be compatible with galactic center excess. To generate a 
singlino mass around 40 to 50 GeV, we vary k = 0.06 ... 0.18, following Eq. (§. 


As mentioned in the previous section, gaugino mass unification correlates the bino, wino, and gluino masses. 
Direct gluino searches set a lower limit of mi /2 > 500 GeV [57]. This leads to a heavy wino mass, out of reach 
for LEP2, and defines a lower bound Mi > 200 GeV. Both, the bino and wino components of the lightest 
neutralino become negligible. To compensate for the missing bino component, the Higgsino component needs 
to be slightly enhanced with respect to the TeV-scale model, leading to the slightly reduced values of /i quoted 
above. 

In the TeV-scale case we fix Ax using Eq. In the NUH-NMSSM this is no longer possible, as Ax is now 
defined at the averaged squark mass scale, where also the Higgs masses are computed. We can estimate that 
for /i = 200 GeV the value of Ax at 10 TeV has to be approximately 8 TeV. Neglecting all contributions but 
Ax itself in Eq.(30), the value of Ax increases to around 8.6 TeV when evaluated at the GUT scale. Erom the 
global analysis we know that A^ tends to have the same sign as Ax. In this case, A^ further increases the 
preferred value of Ax at the GUT scale to around 9 TeV. This fixed value of Ax now translates into preferred 
ranges of A^ and A via the singlet scalar and pseudoscalar squared mass terms, which need to be larger than 
zero. Ghoosing A = 0.25 ... 0.45 gives = (1.5 ... 5) TeV. 


In Eig. we show the results of our SEitter analysis. Just as in the TeV-scale study, we require all 
points to be compatible with direct detection limits, Higgs mass measurements, and the relic density within 
the theoretical uncertainty, as given in Tab. |l| The chargino masses have to be larger than 103 GeV and the 
additional contribution to the Z-width smaller than 2 MeV. All displayed points are consistent with an invisible 
branching ratio of at least 10%. 

As mentioned before, we now study singlino LSPs with a small Higgsino component. On the left hand side of 
Eig. we show the projection onto the k — /i plane. This determines the mass and the composition of the LSP. 
Eor the singlino-like LSP the mass increases with /i and k. The allowed region for larger k corresponds to an LSP 
mass of 50 to 52 GeV while the strip at k ^ O.II corresponds to a neutralino mass of 30 to 40 GeV. In between 
the two regions, the annihilation via the Z-pole becomes too efficient. Eor larger masses the annihilation is too 
weak to predict the measured relic density. In contrast, towards smaller masses a combination of the Ai- and 
Z-channels gives the correct relic density as well as an annihilation cross section compatible with the galactic 
center excess. 

Apart from its mass, the composition of the LSP plays a key role. Eor small k O.I the sum of the Higgsino 
components decreases with increasing /i, starting from 5% at /i = 205 GeV and reaching 20% for /i = 160 GeV. 
This increased active Higgsino component implies a larger coupling to the Higgs, which leads to an increase of 
the invisible branching ratio: for fi < 160 GeV it can reach up to 80%, while for /r > 200 GeV it drops below 
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Figure 6: LSP mass, pseudoscalar mediator mass, and invisible Higgs branching ratio in the NUH-NMSSM. As before, 
we fix tan /3 = 40. Like in Fig. [^all displayed points are compatible with the relic density. Xenon, a chargino mass above 
103 GeV, and the correct SM-like Higgs mass. Moreover, they always have an invisible branching ratio BR(iLi 25 ^ 
Xx) > 10%. The black points are consistent with the galactic center excess. 


the required 10%. However, from the large coupling to the Z there follows a negligible annihilation via the 
pseudoscalar mediator, rendering this region in-compatible with the galactic center excess. Moreover, direct 
detection limits become relevant for a large Z-coupling and exclude points with smaller values of fi. 

In the center panel, we show the A,^ — \ plane for a reduced range of A = 0.30 ... 0.36. This illustrates the 
dependence of the lightest pseudoscalar mass on and A. From Eq.Q now directly follows that the Ai-mass 
increases with A, while it decreases with Once A,^ becomes too large, the pseudoscalar mass squared crosses 
zero, limiting the allowed region. For the galactic center excess av only reaches sufficiently high values around 
the on-shell condition for the Ai funnel. From the discussion of the k — ji plane we know that the mass range 
for neutralinos compatible with the galactic center excess is restricted to = 30 ... 48 GeV. This translated 
into pseudoscalar masses of 60 to 100 GeV. 

In the right panel of Fig. [^the projection on the k — ji plane shows the resulting branching ratio for invisible 
Higgs decays. For the region around k ~ 0.15, A ranges from 0.25 to 0.3, while for ^ 0.11 the allowed range 

for A increases up to 0.45 for /i = 160GeV. Small values of A result in a small Higgsino component, leading 
to an invisible branching ratio of 10 to 30% in the region with ~ 50GeV. For the narrow region the lower 
limit oi 11 = 155GeV in combination with large values of A allow for large invisible branching ratios up to 80%. 

When we consider only points compatible with the GCE we find 


BR(iLi 25 ^ xx) < 40% for ji = 160...200 GeV, (32) 

The maximal found branching ratio of 40% is comparable to the results for /i = 220 GeV in the TeV 
scale NMSSM where the bino component enhances the coupling. Even though the NUH-NMSSM pushes the 
neutralino content to a pure singlino-Higgsino state, we can still find regions that are compatible with the relic 
density, the GCE and a strongly enhanced invisible branching ratio. 


IV. OUTLOOK 

A natural explanation of the Eermi galactic center excess is a light, weakly interacting dark matter particle 
decaying to a pair of bottom quarks through an 5-channel pseudoscalar. The NMSSM is one of the few models 
which predict precisely this process. 

In the NMSSM framework, the galactic center excess as well as the currently observed relic density can be 
accommodated with the help of Z-funnel and Ai-funnel annihilation [26]. Different preferred parameter spaces 
can be linked to LHC searches for trileptons [28] or exotic Higgs searches [30] . We show that for a mixed bino- 
singlino-Higgsino LSP the explanation of the galactic center excess typically predicts large invisible branching 
ratios of the SM-like Higgs boson. In particular for small ji values the invisible branching ratios can reach 50%, 
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testable with Run I Higgs data. Future LHC analyses, sensitive to invisible branching ratios around 3% [2], 
cover a large fraction of Higgs decays to a pair of Hooperons. The preferred NMSSM parameters at the TeV 
scale can also be realized in a unified version of the NMSSM, albeit with larger values of Mi and slightly 
reduced ja. 
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Figure 7: Branching ratio for Hi XX- From left to right we consecutively apply the constraints: (1) mHi 25 — 
(122 ... 128) GeV and BRmv > 10%, (2) m^± > 103 GeV, and (3) = (0.107 ... 0.131). The projections are profiled 

over M 2 , showing the maximal branching ratio. 


Appendix A: Invisible Higgs decays in the MSSM 


For the sake of completion we briefly review the constraints on invisible Higgs decays in the MSSM with a 
SM-like light Higgs. The LSP can be a combination of bino, wino, and Higgsino. First, invisible Higgs decays 
H XX require the LSP to be lighter than 63 GeV. Therefore, at least one of the mass parameters /i, Mi, and 
M 2 has to be around 100 GeV or below. The second ingredient is the size of the coupling. Its form given in 
Eq.(21) requires a mixed Higgsino-gaugino states, which means we again expect to need small values of |/i|. 


The left panel of Fig. [^illustrates the dependence of the invisible branching ratio on fi and Mi. We vary the 
MSSM parameters /i,Mi, and M 2 . For this example we set the sfermion and gluino masses to 2 TeV, so that 
they decouple from the electroweak sector, and as before set tan^^ = 40. The Higgs-sector parameters M^, At 
and now have to be carefully chosen to reproduce the observed Higgs mass. For all points shown we require 
a Higgs invisible branching ratio to be at least 10%. 

Without constraints on the chargino mass and dark matter properties the maximal branching ratio exceeds 
80%. The relevant parameter space is located around |/i| =80 GeV and Mi = 100 ... 150 GeV. Away from this 
region, the LSP is no longer well tempered, reducing the coupling and leading to an invisible branching ratio 
below 10%. 


In the center panel of Fig. [^ we add the LEP limits |22] on the chargino mass. The minimal value of 103 GeV 
translates into a lower bound /i, M 2 > 100 GeV. The lower bound on /i can be seen directly in Eig. The 
constraint on M 2 works indirectly: it excludes wino LSPs lighter than 63 GeV, which means that a light LSP 
has to be mainly bino. This is visible as an upper bound Mi < 80 GeV. 

In addition to the LSP mass, we also have to adjust the couplings. The required bino-Higgsino mixing 
restricts the allowed parameter region to Mi < 80 GeV and |/i| < 200 GeV. Constant invisible branching ratios 
correspond the two half-circles in the center panel of Eig.[^ Without the mass constraints on the chargino and 
the Higgs we would see two approximately circular shapes centered around Mi « 50 GeV and slightly bigger 
values of |/i|. This reflects the preference for a light, well-tempered LSP with roughly equal bino and Higgsino 
fractions. In particular the chargino mass limit simply removes the region with |/i| ^ 100 GeV. The maximum 
invisible branching ratio in the MSSM is 45%. 

In the right panel of Eig. we add the Planck measurement of the LSP relic density. In this configuration 
the annihilation proceeds via an s-channel Z-boson. Planck excludes neutralino masses around 45 GeV, where 
resonant annihilation leads to a too small relic density. Three distinct strips remain compatible with the 
measured relic density: two with neutralino masses between 35 GeV and 40 GeV, and one between 50 GeV 
and 55 GeV. However, the latter is excluded by the XenonlOO. This additional constraint further reduces the 
maximal invisible branching ratio 


BR(Mi25 ^ xx) < 50% 


for 11 = 100 GeV, Mi = 45 GeV , 


(Al) 
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As discussed in Sec. |II A[ in the absence of a light pseudoscalar mediator we do not consider constraints from 
the galactic center excess in the MSSM. 


Finally, we consider Z-decays to neutralinos for neutralino masses smaller than 45 GeV. The corresponding 
partial with adds to the width from Z-decays into neutrinos, whose SM prediction already exceeds the measured 
value of r(Z ^ inv) by 1.9 MeV. While an invisible Higgs branching ratio of 10% only adds an additional 
0.2 MeV to the Z-width, a Higgs branching ratio of 40% can imply an additional Z-decay width of 3 MeV. This 
increases the already existing tension between theory prediction and experimental results for invisible Z-decays. 


[1] O. J. P. Eboli and D. Zeppenfeld, Phys. Lett. B 495, 147 (2000); R. M. Godbole, M. Guchait, K. Mazumdar, 
S. Moretti and D. P. Roy, Phys. Lett. B 571, 184 (2003). 

[2] C. Bernaciak, T. Plehn, P. Schichtel and J. Tattersall, Phys. Rev. D 91, no. 3, 035024 (2015). 

[3] G.-C. Cho, K. Hagiwara, J. Kanzaki, T. Plehn, D. Rainwater and T. Stelzer, Phys. Rev. D 73, 054002 (2006). 

[4] L. Goodenough and D. Hooper, arXiv:0910.2998 [hep-ph]; D. Hooper and L. Goodenough, Phys. Lett. B 697, 412 
(2011); D. Hooper and T. Linden, Phys. Rev. D 84, 123005 (2011); K. N. Abazajian and M. Kaplinghat, Phys. 
Rev. D 86 , 083511 (2012) [Phys. Rev. D 87, 129902 (2013)]; C. Gordon and O. Macias, Phys. Rev. D 88 , no. 
8, 083521 (2013) [Phys. Rev. D 89, no. 4, 049901 (2014)]; O. Macias and C. Gordon, Phys. Rev. D 89, no. 6, 
063515 (2014); K. N. Abazajian, N. Canac, S. Horiuchi and M. Kaplinghat, Phys. Rev. D 90, no. 2, 023526 (2014); 

B. Zhou, Y. F. Liang, X. Huang, X. Li, Y. Z. Fan, L. Feng and J. Chang, arXiv: 1406.6948 [astro-ph.HE]; S. Murgia 
(2014), Talk given on Fifth Fermi Symposium, Nagoya, 20-24 October; B. Bertoni, D. Hooper and T. Linden, 
arXiv: 1504.02087 [astro-ph.HE]. 

[5] P. A. R. Ade et al [Planck Collaboration], arXiv: 1502.01589 [astro-ph.CO]. 

[6] M. Ackermann et al. [Fermi-LAT Collaboration], Phys. Rev. D 89, 042001 (2014). 

[7] [ATLAS Collaboration], Phys. Rev. Lett. 112 , 201802 (2014), etc. 

[8] [CMS Collaboration], CMS-PAS-HIG-13-013; CMS-PAS-HIG-13-018. CMS-PAS-HIG-13-028; arXiv:1310.1002 [hep- 
ex]; CMS-PAS-HIG-14-038; Eur. Phys. J. C 74, 2980 (2014); [CMS Collaboration], arXiv: 1507.00359 [hep-ex]. 

[9] J. R. Ellis, J. F. Gunion, H. E. Haber, L. Roszkowski and F. Zwirner, Phys. Rev. D 39, 844 (1989); M. Drees, 
Int. J. Mod. Phys. A 4, 3635 (1989); J. R. Espinosa and M. Quiros, Phys. Lett. B 279, 92 (1992). U. Ellwanger, 
M. Rausch de Traubenberg and C. A. Savoy, Phys. Lett. B 315, 331 (1993). 

[10] for recent reviews see e.g. S. P. Martin, Adv. Ser. Direct. High Energy Phys. 21, 1 (2010); U. Ellwanger, C. Hugonie 
and A. M. Teixeira, Phys. Rept. 496, 1 (2010); M. Maniatis, Int. J. Mod. Phys. A 25, 3505 (2010). 

[11] for the only other paper with ‘Hooperon’ in the title see A. Alves, S. Profumo, F. S. Queiroz and W. Shepherd, 
Phys. Rev. D 90, no. 11, 115003 (2014). 

[12] T. Daylan, D. P. Finkbeiner, D. Hooper, T. Linden, S. K. N. Portillo, N. L. Rodd and T. R. Slatyer, arXiv: 1402.6703 
[astro-ph.HE]; F. Calore, 1. Cholis and C. Weniger, JCAP 1503, 038 (2015); F. Calore, 1. Cholis, C. McCabe and 

C. Weniger, Phys. Rev. D 91, no. 6, 063003 (2015). 

[13] H. Goldberg, Phys. Rev. Lett. 50, 1419 (1983); M. Drees and M. M. Nojiri, Phys. Rev. D 47, 376 (1993); J. R. Ellis, 
J. S. Hagelin, D. V. Nanopoulos, K. A. Olive and M. Srednicki, Nucl. Phys. B 238, 453 (1984); G. Jungman, 
M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 (1996); G. Bertone, D. Hooper and J. Silk, Phys. Rept. 
405, 279 (2005). 

[14] for an exceptionally instructive and beautiful set of plots on relic neutralinos see J. Bramante, P. J. Fox, A. Martin, 

B. Ostdiek, T. Plehn, T. Schell and M. Takeuchi, Phys. Rev. D 91, no. 5, 054015 (2015). 

[15] P. Agrawal, B. Batch, P. J. Fox and R. Harnik, JCAP 1505, no. 05, 011 (2015); A. Achterberg, S. Caron, L. Hendriks, 
R. Ruiz de Austri and C. Weniger, arXiv: 1502.05703 [hep-ph]. 

[16] S. Henrot-Versille, R. Lafaye, T. Plehn, M. Rauch, D. Zerwas, S. Plaszczynski, B. Rouille d’Orfeuil and M. Spinelli, 
Phys. Rev. D 89, no. 5, 055017 (2014). 

[17] T. Cohen and J. G. Wacker, JHEP 1309, 061 (2013). 

[18] see e.g. A. Fowlie, K. Kowalska, L. Roszkowski, E. M. Sessolo and Y. -L. S. Tsai, Phys. Rev. D 88, 055012 (2013); 

C. Strege, G. Bertone, F. Feroz, M. Fornasa, R. Ruiz de Austri and R. Trotta, JCAP 1304, 013 (2013); B. Bhat- 
tacherjee, M. Chakraborti, A. Chakraborty, U. Chattopadhyay, D. Das and D. K. Ghosh, arXiv: 1305.4020 [hep-ph]; 
M. Cahill-Rowley, J. Hewett, A. Ismail and T. Rizzo, arXiv:1308.0297 [hep-ph]; O. Buchmiiller, R. Cavanaugh, 
M. Citron, A. De Roeck, M. J. Dolan, J. R. Ellis, H. Flacher and S. Heinemeyer et al, Eur. Phys. J. C 72, 2243 
(2012); etc. 

[19] K. Griest and D. Seckel, Phys. Rev. D 43, 3191 (1991); S. Mizuta and M. Yamaguchi, Phys. Lett. B 298, 120 
(1993); J. R. Ellis, T. Falk and K. A. Olive, Phys. Lett. B 444, 367 (1998). 

[20] C. Boehm, A. Djouadi and M. Drees, Phys. Rev. D 62, 035012 (2000). 






23 


[21] P. Binetmy, G. Girardi and P. Salati, NucL Phys. B 237, 285 (1984); S. Mizuta and M. Yamaguchi, Phys. Lett. B 
298, 120 (1993); 

[22] ALEPH Gollaboration, Phys. Lett. B526 (2002) 206, Phys. Lett. B583 (2004) 247; DELPHI Gollaboration, Eur. 
Phys. J. G31 (2003) 421-479; L3 Gollaboration, Phys. lett B580 (2004) 37; OPAL Gollaboration, Eur. Phys. J 
G32 (2004) 453-473; ALEPH Gollaboration, Phys. Lett. B537 (2002) 5-20; OPAL Gollaboration, Phys. Lett. B545 
(2002) 272-284; ALEPH Gollaboration, Phys.Lett. B533 (2002) 223-236; L3 Gollaboration, Phys. Letter B471 (999) 
308-320. LEP2 SUSY Working Group, ALEPH, DELPHI, L3 and OPAL experiments, http://lepsusy.web.cern. 
ch/lepsusy 

[23] D. K. Ghosh, S. Mondal and 1. Saha, JGAP 1502, no. 02, 035 (2015); A. Berlin, S. Gori, T. Lin and L. T. Wang, 
arXiv:1502.06000 [hep-ph]. 

[24] A. Berlin, D. Hooper and S. D. McDermott, Phys. Rev. D 89, no. 11, 115022 (2014); M. Abdullah, A. DiEranzo, 
A. Rajaraman, T. M. P. Tait, P. Tanedo and A. M. Wijangco, Phys. Rev. D 90, no. 3, 035004 (2014); G. Balazs, 
T. Li, G. Savage and M. White, arXiv: 1505.06758 [hep-ph]. 

[25] for general aspects of the singlet/singlino in the NMSSM see e.g. U. Ellwanger and G. Hugonie, Eur. Phys. J. G 
5, 723 (1998); A. Dedes, G. Hugonie, S. Moretti and K. Tamvakis, Phys. Rev. D 63, 055009 (2001); S. Y. Ghoi, 
D. J. Miller and P. M. Zerwas, Nucl. Phys. B 711, 83 (2005); V. Barger, P. Langacker and G. Shaughnessy, Phys. 
Lett. B 644, 361 (2007); K. Gheung and T. J. Hon, Phys. Lett. B 674, 54 (2009); O. Stal and G. Weiglein, JHEP 
1201, 071 (2012) D. G. Gerdeno, P. Ghosh, G. B. Park and M. Peiro, JHEP 1402, 048 (2014); U. Ellwanger and 
A. M. Teixeira, JHEP 1410, 113 (2014). 

[26] M. Gahill-Rowley, J. Gainer, J. Hewett and T. Rizzo, JHEP 1502, 057 (2015); T. Gherghetta, B. von Harling, 
A. D. Medina, M. A. Schmidt and T. Trott, Phys. Rev. D 91, 105004 (2015); R. Enberg, S. Munir, G. P. d. 1. Heros 
and D. Werder, arXiv: 1506.05714 [hep-ph]; J. Gao, L. Shang, P. Wu, J. M. Yang and Y. Zhang, arXiv: 1506.06471 
[hep-ph], 

[27] J. Huang, T. Liu, L. T. Wang and F. Yu, Phys. Rev. D 90, no. 11, 115006 (2014). 

[28] J. Guo, J. Li, T. Li and A. G. Williams, Phys. Rev. D 91, no. 9, 095003 (2015); J. Gao, L. Shang, P. Wu, J. M. Yang 
and Y. Zhang, Phys. Rev. D 91, no. 5, 055005 (2015). 

[29] X. J. Bi, L. Bian, W. Huang, J. Shu and P. E. Yin, arXiv: 1503.03749 [hep-ph]. 

[30] G. Gheung, M. Papucci, D. Sanford, N. R. Shah and K. M. Zurek, Phys. Rev. D 90, no. 7, 075011 (2014). 

[31] J. Kozaczuk and T. A. W. Martin, JHEP 1504, 046 (2015). 

[32] for a Higgs portal linked to invisible decays see e.g. V. Silveira and A. Zee, Phys. Lett. B 161, 136 (1985); J. Mc¬ 
Donald, Phys. Rev. D 50, 3637 (1994); G. P. Burgess, M. Pospelov and T. ter Veldhuis, Nucl. Phys. B 619, 709 
(2001); R. Schabinger and J. D. Wells, Phys. Rev. D 72, 093007 (2005) B. Patt and E. Wilczek, hep-ph/0605188. 
R. Barbieri, T. Gregoire and L. J. Hall, hep-ph/0509242; M. H. G. Tytgat, PoSIDM 2010, 126 (2011); G. Englert, 
T. Plehn, D. Zerwas and P. M. Zerwas, Phys. Lett. B 703, 298 (2011); M. Pospelov and A. Ritz, Phys. Rev. D 
84, 113001 (2011) X. -G. He and J. Tandean, Phys. Rev. D 84, 075018 (2011); P. J. Eox, R. Harnik, J. Kopp 
and Y. Tsai, Phys. Rev. D 85, 056011 (2012); 1. Low, P. Schwaller, G. Shaughnessy and G. E. M. Wagner, Phys. 
Rev. D 85, 015009 (2012); G. Englert, T. Plehn, M. Rauch, D. Zerwas and P. M. Zerwas, Phys. Lett. B 707, 512 
(2012); A. Djouadi, O. Lebedev, Y. Mambrini and J. Quevillon, Phys. Lett. B 709, 65 (2012); B. Batch, S. Gori 
and L. T. Wang, JHEP 1206, 172 (2012); S. Back, P. Ko and W. 1. Park, Phys. Rev. D 90, no. 5, 055014 (2014). 

[33] J. R. Ellis, J. E. Gunion, H. E. Haber, L. Roszkowski and E. Zwirner, Phys. Rev. D 39, 844 (1989); M. Drees, 
Int. J. Mod. Phys. A 4, 3635 (1989); E. Eranke and H. Eraas, Int. J. Mod. Phys. A 12 , 479 (1997); D. J. Miller, 
R. Nevzorov and P. M. Zerwas, Nucl. Phys. B 681, 3 (2004); U. Ellwanger and G. Hugonie, Phys. Lett. B 623, 93 
(2005); G. Degrassi and P. Slavich, Nucl. Phys. B 825, 119 (2010). 

[34] M. Masip, R. Munoz-Tapia and A. Pomarol, Phys. Rev. D 57, R5340 (1998); M. A. Diaz, T. A. ter Veldhuis 
and T. J. Weiler, Phys. Rev. D 54, 5855 (1996); U. Ellwanger and G. Hugonie, Eur. Phys. J. G 25, 297 (2002); 
M. D. Goodsell, K. Nickel and E. Staub, Phys. Rev. D 91, no. 3, 035021 (2015). 

[35] [ATLAS and GMS Gollaborations], Phys. Rev. Lett. 114, 191803 (2015). 

[36] for a short history of MSSM studies see e.g. M. S. Garena, J. R. Espinosa, M. Quiros and G. E. M. Wagner, Phys. 
Lett. B 355, 209 (1995); H. E. Haber, R. Hemphing and A. H. Hoang, Z. Phys. G 75, 539 (1997); S. Heinemeyer, 
W. Hollik and G. Weiglein, Eur. Phys. J. G 9, 343 (1999); G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich and 
G. Weiglein, Eur. Phys. J. G 28, 133 (2003); P. Kant, R. V. Harlander, L. Mihaila and M. Steinhauser, JHEP 1008, 
104 (2010); T. Hahn, S. Heinemeyer, W. Hollik, H. Rzehak and G. Weiglein, Phys. Rev. Lett. 112 , no. 14, 141801 
(2014); E. Bagnaschi, R. V. Harlander, S. Liebler, H. Mantler, P. Slavich and A. Vicini, JHEP 1406, 167 (2014); 
G. Degrassi, S. Di Vita and P. Slavich, Eur. Phys. J. G 75, no. 2, 61 (2015). 

[37] G.W. Bennett et ah, Phys. Rev. Lett. 89, 101804 (2002); Erratum ibid. Phys. Rev. Lett. 89, 129903 (2002); G.W. 
Bennett et ah, Phys. Rev. Lett. 92, 161802 (2004); G.W. Bennett et ah, Phys. Rev. D73, 072003 (2006); Particle 
Data Group http: //pdg. Ibl .gov/2014/reviews/rpp2014-rev-g-2-muon-anom-mag-moment .pdf 

[38] GLEG Gollaboration, Phys. Rev. Lett. 87, 251807 (2001); BaBar Gollaboration, Phys. Rev. D86, 112008 (2012); 
Belle Gollaboration, Phys. Rev. Lett. 103, 241801 (2009); Particle Data Group http://pdg.lbl.gov/2014/reviews/ 
rpp2014-rev-b-meson-prod-decay.pdf 

[39] [LHGb Gollaboration], Phys. Rev. Lett. 110, 021801 (2013); GMS Gollaboration, Phys. Rev. Lett. Ill, 101804 


24 


(2013); CMS and LHCb Collaborations, Nature 522, 68 (2015). 

[40] D. Asner et al. [Heavy Flavor Averaging Group Collaboration], arXiv:1010.1589 [hep-ex]. 

[41] ALEPH Collaboration, EPJ C14 1; Delphi Collaboration, EPJ C16 371; L3 Collaboration, EPJ C16 1; OPAL 
Collaboration, EPJ C19 587; Particle Data Group: http://pdg.lbl.gov/2014/listings/rpp2014-list-z-boson. 
pdf 

[42] Tevatron Electroweak Working Group [CDE and DO Collaborations], arXiv: 1407.2682 [hep-ex]; [ATLAS and 
CDE and CMS and DO Collaborations], arXiv: 1403.4427 [hep-ex]; K. Kroninger, A. B. Meyer and P. Uwer, 
arXiv: 1506.02800 [hep-ex]; Particle Data Group: http://pdg.lbl.gov/2014/reviews/rpp2014-rev-top-quark. 
pdf 

[43] E. Aprile et al [XENONIOO Collaboration], Phys. Rev. Lett. 107, 131302 (2011); E. Aprile et al [XENONIOO 
Collaboration], Phys. Rev. Lett. 109, 181301 (2012). We use the conservative XenonlOO results, but in our bottom 
line the only effects from including the Lux limits is a shift in /x and its correlations. 

[44] R. Lafaye, T. Plehn and D. Zerwas, arXiv:hep-ph/0404282; R. Lafaye, T. Plehn, M. Rauch and D. Zerwas, Eur. 
Phys. J. C 54, 617 (2008). 

[45] C. Adam, J. -L. Kneur, R. Lafaye, T. Plehn, M. Rauch and D. Zerwas, Eur. Phys. J. C 71, 1520 (2011); E. Turlay, 
R. Lafaye, T. Plehn, M. Rauch and D. Zerwas, J. Phys. G 38, 035003 (2011). 

[46] for the similar Pitting approach see e.g. P. Bechtle, K. Desch and P. Wienemann, arXiv:hep-ph/0412012; P. Bechtle, 
T. Bringmann, K. Desch, H. Dreiner, M. Hamer, C. Hensel, M. Kramer and N. Nguyen et al, JHEP 1206 (2012) 
098. 

[47] A. Hocker, H. Lacker, S. Laplace and E. Le Diberder, Eur. Phys. J. C 21, 225 (2001). 

[48] U. Ellwanger, J. E. Gunion and C. Hugonie, JHEP 0502 (2005) 066 [hep-ph/0406215]. U. Ellwanger and C. Hugonie, 
Comput. Phys. Commun. 175, 290 (2006); G. Belanger, E. Boudjema, C. Hugonie, A. Pukhov and A. Semenov, 
JCAP 0509, 001 (2005); U. Ellwanger and C. Hugonie, Comput. Phys. Commun. 177, 399 (2007). 

[49] G. Belanger, E. Boudjema, P. Brun, A. Pukhov, S. Rosier-Lees, P. Salati and A. Semenov, Comput. Phys. Commun. 
182, 842 (2011). 

[50] R. Lafaye, T. Plehn, M. Rauch, D. Zerwas and M. Diihrssen, JHEP 0908, 009 (2009); M. Klute, R. Lafaye, T. Plehn, 
M. Rauch and D. Zerwas, Phys. Rev. Lett. 109, 101801 (2012); D. Lopez-Val, T. Plehn and M. Rauch, JHEP 1310, 
134 (2013). 

[51] A. Djouadi, J. Kalinowski and M. Spira, Comput. Phys. Commun. 108, 56 (1998); M. Miihlleitner, A. Djouadi and 
Y. Mambrini, arXiv:hep-ph/0311167; A. Djouadi, M. M. Miihlleitner and M. Spira, Acta Phys. Polon. B 38, 635 
(2007). 

[52] A. Djouadi, J. -L. Kneur and G. Moultaka, Comput. Phys. Commun. 176, 426 (2007); G. Brooijmans et al, 
arXiv:1203.1488 [hep-ph]. 

[53] S. Heinemeyer, W. Hollik, A. M. Weber and G. Weiglein, JHEP 0804, 039 (2008). A. Weber, private communication. 

[54] T. Corbett, O. J. P. Eboli, D. Goncalves, J. Gonzalez-Eraile, T. Plehn and M. Rauch, arXiv:1505.05516 [hep-ph]. 

[55] L. E. Ibanez, Phys. Lett. B 118, 73 (1982); J. R. Ellis, D. V. Nanopoulos and K. Tamvakis, Phys. Lett. B 121, 
123 (1983); L. Alvarez-Gaume, J. Polchinski and M. B. Wise, Nucl. Phys. B 221, 495 (1983); K. Inoue, A. Kakuto, 
H. Komatsu and S. Takeshita, Prog. Theor. Phys. 68 , 927 (1982) [Erratum-ibid. 70, 330 (1983)]; A. H. Chamseddine, 
R. Arnowitt and P. Nath, Phys. Rev. Lett. 49, 970 (1982). 

[56] A. Djouadi, U. Ellwanger and A. M. Teixeira, Phys. Rev. Lett. 101, 101802 (2008). 

[57] [ATLAS Collaboration], JHEP 1409, 176 (2014). 


